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Methods of Predicting the Atmospheric 
Bending oi Radio Rays 1 
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Three methods for predicting the bending of radio rays when the refractive index profile 
above the surface layer is unknown have been developed recently by the authors. These 
methods are: a statistical technique for refraction at high initial elevation angles, estimation 
of bending from an exponential model of atmospheric refractive index, and a modification 
of the exponential model to account for the heavily weighted effects of anomalous initial 
refractive index gradients at small initial elevation angles. Each model is dependent upon 
the value of the refractive index at ground level or, in the case of superrefraction, the addi- 
tional knowledge of the refractive index gradient next to the earth's surface. Each method 
works best in a particular range of initial elevation angles <>r meteorological conditions. 
The heighl and angular ranges of application of each met hod are checked by comparison wit li 
values obtained from 77 diverse refractive index profiles representative of wide climatic 

variation. It is found that the use of the best of the three methods will always result in a 
prediction of the total atmospheric bending within 10 percent for initial elevation angles 
from zero to 10 milliradians and to within I percent for initial elevation angles greater than 
1 7 milliradians (~ 1 deg). 
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Glossary of Terms 

n=\ he radio refractive index. 

n s =the value of n at the earth's surface. 

i\T=the radio refractivity, N= (n— 1)X10'\ 
JV s =value of A r at the earth's surface. 
AN= difference between N s and the N value 
at one km above the surface, —AN= 

(dN/dh)o=gT&dient of N with respect to 
dN/dh, evaluated at the 
surface. 

r=vector radius from the center 
earth. 

r =radial distance from the center 
earth's surface. 

h = height above the surface, h=r—r . 
#— elevation angle of a radio ray, 

(acute) angle between the tangent to 
the ray path and the local horizontal 
(i.e., perpendicular to the radius 
vector.) 

O = the value of 6 at the ray path origin 
(transmitting or receiving point). 

6 p =ihe angle of penetration for a radio 
duct, i.e., the smallest value of O for 
which the radio ray will not be 
trapped, or conversely, the largest O 
for which the ray will be trapped. 
r = the angular refraction, or bending, of a 

radio ray. 
e — elevation angle error, the angular differ- 
ence between O and the true elevation 
angle to a target at a given point on 
the ray path. 
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i Contribution from Central Radio Propagation Laboratory, National Bureau 
of Standards, Boulder, Colo. 



1. Introduction 

Recent years have seen considerable activity in 
the evaluation of refraction effects in the troposphere. 
Schulkin [l] 2 outlined a simple method for refraction 
calculations and applied it to determine the mean 
refraction expected in arctic, temperate, and tropical 
climates. Fannin and Jehu [2] made an extensive 
analysis of elevation angle errors expected in various 
air masses and geographic locations for initial eleva- 
tion angles in excess of 3 (\^j: above the horizontal. 
In a series of papers [3,4,5,6] the present authors 
have examined atmospheric refractive index struc- 
ture, the effect of this structure upon radio-ray 
refraction and have evolved methods of estimating 
the refraction of radio rays for all initial elevation 
angles. These methods are unique in that the\ 
depend only upon a knowledge of the refractive 
index at the earth's surface or, in the case of super- 
refraction, the gradient of the refractive index in the 
earth-boundary layer. Thus distributions of eleva- 
tion angle error, angular bending, and other refrac- 
tion effects may be determined for the majority of 
practical applications by simple reference to distri- 
butions of the surface value of the refractive index 
such as those foi the United States [7], It is assumed 
of course, that these various methods will be applied 
only when either details of the actual refractive 
index profile are unknown or it is impractical to 
obtain these details. It is under this assumption 
that this paper tests the relative accuracy of these 
various prediction methods to arrive at a delineation 
of conditions under which each method works best. 



1 Figures in brackets indicate the literature references at the end of this paper. 
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2. Theory and Background 

A basic measure of atmospheric refraction effects 
on high frequency radio propagation is the bending, 
or angular refraction, of individual rays. By as- 
suming the refractive index to be a function only of 
height above the surface of a smooth, spherical earth 
rays can be traced using SnelPs law in the following 
form [1]: 



nr cos 6=n s r cos 6 Q 



(i) 



The geometry is shown in figure 1. The equation 
for the ray bending, r, can be obtained from eq (1) 
as [1, 9] 

'... = - J^cot* 



where 



cos 



Jn s n 
_i f Vo cos d \ 



and cot 6 is given to a very high degree of approxima- 
tion by: 



cos i 



: V sii 



cot 0= A / S m 2 o +— -2(Af s -AOX10- 6 cos 2 6 ' 



It can be seen from inspection of eq (2) and (3) 
that in order to evaluate t directly it is necessary to 
know the refractive index of the atmosphere, n as a 
function of height, h. 

How one evaluates eqs (2) and (3) depends upon 
the availability of data and computing facilities. If 
details of the refractive index profile are available 




Figure 1. Geometry of radio-ray refraction. 



from either refractometer or radiosonde ascents, then 
one could evaluate refraction effects with the aid of 
of an electronic computer or simplified graphical 
techniques [9]. However, frequently the desired 
refractive index profile data are simply not available 
and one must fall back upon less exact methods of 
estimating r. One may use average values of r such 
as given by Schulkin or the distributions of refraction 
effects given by Fannin and Jehn. However, simple, 
standard measurements of pressure, temperature, 
and humidity at the earth's surface are almost 
always available and may be used to estimate tropo- 
spheric refraction effects. In the sections that follow, 
the evaluation of the integral for r will be examined 
from three different viewpoints: 

1. Simplification of eqs (2) and (3), to permit 
evaluation of r without knowledge of details of the 
actual refractive index profile. l 

2. Evaluation of r for actual observed refractive 
index profiles and the statistical reduction of the data 
so derived into a function of some observable pa- 
rameter (e.g., N s ). 

3. Construction of an analytic model of n(r) for 
normal conditions, thus yielding expected values of 
r by direct integration of the model. 

3. r as a Function of N s 

As a first approximation towards evaluation of r 
without detailed knowledge of n(r), consider the 
integration by parts of eq (2): 



T s,l z 



or: 



pick 

Jn n 



]n 1 

/'cot 2 

+ In {n}d (cot 6) 

J cot n 



t s i=ln {n s } cot0 o — In {^Jcot 0i~ In {n} csc 2 6dd 

Now since n^l+A^XlO" 6 and 2VX1CT 6 <5X10- 4 
then 



In [n] =NX 10- 6 { 1 - % (NX 10" 



.}^7VX10"^ 



with this approximation (i.e., a maximum error of 
less than 0.2%) the above equation becomes 



T s ,l< 



^iV s XlO- 6 cot0 o 



-r 



-iVjXlO^cotfl! 

A/X10- 6 csc 2 ddd for O ±O. 



(4) 



The integral in eq (4) has been found to contribute 
no more than 3 percent to the value of r S;CO for an 
initial elevation angle of 10 deg or greater [1] while 
the second term is zero due to N m =0. Thus the 
first term of (4) forms an approximation to r S)CO/ 
which is asymptotic to the true value of r s>00 as O 
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approaches 90 cleg. It may be shown that for nor- 
mal conditions and all heights the integral in (4) is 
essentially independent of N s for 6 o y>l7 mr (~1 
deg); the term N x cot di tends to be constant; thus 
(4) reduces to a linear equation 



Ts^bxNs+ai. 



(5) 



The form of eq (5) is very attractive, since it im- 
plies two things: 

1. t Sj i may be predicted with some accuracy as a 
function only of N s (hi and O constant), a parameter 
which may be observed from simple surface meas- 
urements of the common meteorological elements of 
temperature, pressure, and humidity. 

2. The simple linear form of the equation indicates 
that, given a large mass of observed r S( i versus N s 
for many values of h and 6 0} the expected (or best 
estimate) values of b and a can be obtained by the 
standard method of statistical linear regression. 

The method of attack indicated by implication 
number two has, in fact, been carried out by the 
authors [3]. The results show that for h x — oo the 
method is accurate to within ±3 percent of the true 
value (as an rms error) for initial elevation angles as 
small as 1 deg. The accuracy of this method and 
the following methods, will be examined for 0.1 
km<A< oo in the following sections. 

4. Exponential Model 

The development of a model of N(h) to describe 
the normal behavior of atmospheric N as a function 
of N s and height has received a considerable amount 
of treatment in the past. One of the simplest, and, 
at the same time, most accurate models which has 
emerged from these studies is that in which N(h) has 
an exponential decrease with height [4, 5, 10], 



N(h)=N s exp { — c e h}. 



(6) 



One of the earliest applications of this particular 
model has been attributed by Garfinkel to Sir Isaac 
Newton, who used the exponential form in a study 
of astronomical refraction [11]. 



If it is assumed that N(h) is indeed an exponential 
function of height, then the gradient of N(h) would 
also be an exponential function of height. The most 
extensive source of data with which to evaluate the 
coefficients in the exponential is that of AN (the 
value of N at 1 km minus the surface value, N s ) 
which has received wide application in radio propa- 
gation problems [12]. Thus we would expect 



AN 
Ah 



=ki exp {— k 2 h} 



(7) 



to take the form 



AN=h exp {— k 2 } 

for our special case of Ah=h=l km. Examination 
of the AN data soon revealed that k 2 was dependent 
upon N s i.e., the higher the surface value of N the 
greater the expected drop in iVover 1 km. Exami- 
nation of the data indicated that 

k 2 =k z N s 

and the resultant equation, 

AA^exp {-k 3 N s } (8) 

was solved by least squares. The least squares de- 
termination was facilitated by converting (8) to the 
form 

In {AN^-ksNs+lnh (9) 

or, in words, expressing the natural logarithm of AN 
as a linear function of N s . The values of k x and k 3 
were established from some 888 sets of 8-year means 
of AN and N s from 45 U.S. weather stations. The 
results of this study are shown graphically in figure 
2_where the least squares exponential fit of AA r and 
N s is given by: 



-AN=7.32 exp {0.005577iV s 
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Figure 2. Regression of In |aa| upon Ns. 
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From this equation the CRPL Exponential Reference 
Atmosphere [5] was determined, the profiles being de- 
fined by the following equations: 



N=N S exp {— c e h} 



c e =ln 



\N s +ANf 



(11) 



Ray tracings have been computed for this model 
covering more than the normal range of N s , and the 
results may be used, either in tabular or graphical 
form, to predict r for any normal combination of N S} 
So, and height [5]. 

5. Initial Gradient Correction Method 

The importance of the initial gradient in radio 
propagation where the initial elevation angle of a ray 
path is near zero has long been recognized. For ex- 
ample, if dN/dh=l/r then At= qo, an expression of 
the fact that the ray path will travel at a constant 
height above the earth's surface. This is called 
ducting, or trapping of the radio ray. The effect of 
anomalous initial iV-gradients on ray propagation at 
elevation angles near zero, and for gradients less than 
ducting, (\dN/dh\<157/km, or dn/dh>- 157 /km), 
may also be quite large. A method has been pre- 
sented [4] for correcting the predicted refraction 
(from the exponential reference atmosphere) to ac- 
count for anomalous initial ^/-gradients, assuming 
that the actual value of the initial gradient is known. 
The result is, 

r h =r h (N Sl d ) + [r m (Ntd Q )-r m (N s ,d )} (12) 

where T h (N s ) = T at height h, for the exponential ref- 
erence atmosphere corresponding to N s and N* is 
the N s for the exponential reference atmosphere hav- 
ing the same initial gradient as that observed; r m is 
r at 100 meters height. 

This procedure has the effect of correcting the pre- 
dicted bending by assuming that the observed initial 
gradient exists throughout a surface layer 100 meters 
thick, calculating the bending at the top of the 100- 
meter-thick layer, and assuming that the atmosphere 
behaves according to the exponential reference profile 
corresponding to the observed value of N s for all 
heights above 100 meters. This approach has proven 
quite successful in predicting r for initial elevation 
angles under 10 mr, and will, of course, predict duct- 
ing when it occurs. 

6. Analysis of the Accuracy of the Predic- 
tion Methods 

A test sample of ray bendings for the range of N- 
profiles likely to be encountered was prepared from 
77 refractivity profiles derived from radiosonde ob- 
servations. These 77 profiles represent both normal 
and extreme refractivity profiles for 13 climatically 
diverse locations in the United States and are actu- 
ally representative of a nearly worldwide range of 
conditions [13]. Values of r were calculated by 



numerical integration for values of O from to 900 
mr. (For a thorough discussion of the ray-tracing 
techniques employed see reference [5].) The results 
of this general refraction study provided a large mass 
of data for checking the accuracy of each of the pre- 
diction methods. 

The relative accuracy of the exponential model 
and the initial gradient correction method were tested 
by predicting the bending at particular heights and 
initial elevation angles for the 77 sample profiles and 
finding the rms error of prediction for each case. 
Since 13 of the 77 profiles had surface ducts they 
could be used only for elevation angles greater than 
the angle of penetration for each case, given by: 



d p =-yj0.2[156.9-(dN/dh) ]X10 Q radians (13) 

where (dN/dh) is the observed initial gradient of TV 
per km, assumed to extend over 100 m. Thus there 
were only 64 profiles analyzed at 6 =0, and 77 at O = 
10 mr. These same data were also used to derive 
the regression of r upon N s for various heights and 
initial elevation angles. The scatter of points about 
the regression line was then used as an estimate of the 
minimum rms error that would be expected from any 
of the three prediction methods. Note carefully, 
however, that the regression lines are a "best fit" to 
the test data. Although a future sample of data 
would presumably yield similar rms deviations, the 
possibility exists that the present sample data are 
systematically biased, in this case there would be an 
additional error encountered in practice due to this 
bias. 

The primar3 r purpose of the present analysis is to 
determine over which regions of height and 6 each 
method will give the best results. It was found that 
the statistical correlation method is more accurate 
from 1 deg to vertical incidence and that for all alti- 
tudes in this range of O it is the most accurate of the 
three methods. At 6 smaller than 1 deg, and espe- 
cially at 6 <C10 mr, the methods based on the ex- 
ponential reference atmosphere, particularly the cor- 
rected exponential reference atmosphere, are more 
accurate. Figure 3 illustrates these conclusions by 
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Figure 3. RMS error of predicting r at various heights as a 
percent of mean t, excluding super refractive profiles. 
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comparing the rms error of prediction as a percent of 
mean r for all three methods at three different heights 
and over a large range of O , excluding all superrefrac- 
tive profiles. A superrefractive profile is here defined 
as one with an initial gradient of N in excess of 100 
iV-units per km, i.e., (dN/dh) < — 100/km. 

It is evident from figure 3 that the percentage error 
of predicting r decreases with increasing thickness of 
the atmosphere through which the ray passes. This 
is as one would expect since the value of the refrac- 
tive index becomes less variable with increasing 
height and, one might say, the upper limit of integra- 
tion of (2) becomes more a function of the lower 
limit, N s . Further, the sensitivity of refraction 
effects to low-level profile anomalies for small values 
of O is reflected by the relatively small percentage 
error of the corrected exponential model for flo^lO 
mr. 

It must be remembered that the statistical regres- 
sion technique is only apparently superior to methods 
using the exponential reference atmosphere since, by 
definition, it must have a minimum rms error pro- 
vided that r is a linear function of N s . Although the 
difference in the percentage rms error between these 
two methods appears quite large, the actual rms 
error of the exponential reference atmosphere is 
either less than or within 0.1 mr of that of the statis- 
tical correlation technique for 6 >7 mr, thus, for 
practical considerations, indicating no clear-cut su- 
periority of one technique over the other; this is 
illustrated in figure 4. 
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Figure 4. Difference between rms error of predicting r at 70 
km for the exponential reference atmosphere and the statistical 
correlation method. 

One may evaluate the minimum rms error expected 
in predicting r by these three methods by defining a 
composite prediction method that utilizes the best 
of the three methods in each range of O and height 
to yield a minimum overall error. 

The numerical value of the rms deviations for the 
optimized composite of the three methods is shown 
on figure 5 for heights of 1, 3, and 70 km and all 
profiles. It is seen that the maximum error for any 
case is about 2 mr and decreases to 1 mr or less for 
#o>10 mr for all height increments. The error for 
the total bending case, h=70 km, drops below an 
rms value of 1 mr for 6 >5 mr. 
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Figure 5. RMS error of predicting r, composite prediction 
method. 



7. Conclusions 

The present study appears to indicate that: 

1 . The statistical regression technique is an ade- 
quate solution to the bending problem for all O 
larger than about 10 mr, and all heights from 1 km 
up. 

2. The exponential reference atmosphere is equally 
as good as the statistical technique over the same 
ra nge of t) and height. 

3. The initial gradient correction method is very 
useful for <^1() mr, and can be used to extend 
predictions of r down to a d of for any model which 
vie Ids N as a function of height. 

These conclusions indicate that the reader who 
desires the quick evaluation of some refraction effect 
should consult the rather extensive tables of the 
CRPL Exponential Reference Atmosphere [5]. 
These tables would allow, for instance, the deter- 
mination of the elevation angle error as a function 
of O , N at the earth's surface, and the radar range. 
The reader desiring the convenience of the statistical 
regression technique for estimating either r or the 
elevation angle error may obtain the necessary 
statistical parameters from the literature [14]. 



The authors express their gratitude to Mrs. B. J. 
Weddle, Mrs. M. A. Fischer, and Mrs. G. M. 
Richmond for their aid in the calculations of this 
study. 

8. References 

[1] M. Schulkin, Average radio-ray refraction in the lower 
atmosphere, Proc. IRE 40, 554 (1952) 

[2] B. M. Fannin and K. H. Jehn, A study of radar eleva- 
tion-angle errors due to atmospheric refraction, IRE 
Trans. AP-5, 71 (1957). 

[3] B. R. Bean and B. A. Cahoon, The use of surface weather 
observations to predict the total atmospheric bending 
of radio waves at small angles, Proc. IRE 45, 1545 
(1957). 

[4] B. R. Bean and G. D. Thayer, On models of the atmo- 
spheric radio refractive index, Proc. IRE 47, 740 (1959) 

[5] B. R. Bean and G. D. Thayer, CRPL exponential refer- 
ence atmosphere, NBS Monograph 4 (1959). 



491 



[6] B. R. Bean and B. A. Cahoon, Effect of atmospheric 
horizontal inhomogeneity upon ray tracing, J. Re- 
search NBS 63D, 287 (1959). 

[7] B. R. Bean, J. D. Horn, and A. M. Ozanich, Jr., Climatic 
charts and data of the radio refractive index for the 
United States and the world, NBS Monograph (in 
press) . 

[8] W. M. Smart, Spherical astronomy, chap. 3 (Cambridge 
University Press, London, 1931). 

[9] S. Weisbrod and L. J. Anderson, Simple methods for 
computing tropospheric and ionospheric refractive 
effects on radio waves, Proc. IRE 47, 1770 (1959). 

[10] J. R. Bauer, W. C. Mason, and F. A. Wilson, Radio re- 
fraction in a cool exponential atmosphere, Lincoln 
Laboratory, MIT, Tech. Rpt. 186 (27 Aug. 1958). 



[11] Boris Garfinkel, An investigation in the theory of 
astronomical refraction, The Astronomical Journal 50, 
169 (1944). 

[12] B. R. Bean and F. M. Meaney, Some applications of the 
monthly median refractivitv gradient in tropospheric 
propagation, Proc. IRE 43,"l419 (1955). 

[13] B. R. Bean and G. D. Thayer, Comments on the limits 
to the utilization of the refractive index at ground level 
as a radio-meteorological parameter, To be published 
in Proc. IRE in reply to a letter by P. Misme. 

[14] B. R. Bean, B. A. Cahoon, and G. D. Thayer, Tables 
for the statistical prediction of radio ray bending and 
elevation angle error using the surface values of the 
radio refractive index, NBS Tech. Note 44 (in press). 

(Paper 64D5-84) 



492 



